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We present a search for direct CP violation in B^ — > J/ipK ± (n ) decays. The event sample is se- 
lected from 2.8 ftr 1 of pp collisions recorded by DO experiment in Run II of the Fermilab Tevatron 
Collider. The charge asymmetry A C p(B+ -> J/ipK+) = +0.0075 ± 0.0061(stat.)±0.0027(syst.) 
is obtained using a sample of approximately 40 thousand B ± — > J/tpK ± decays. The achieved 
precision is of the same level as the expected deviation predicted by some extensions of the 
standard model. We also measured the charge asymmetry Acp(B + — > J/ipir + ) = —0.09 ± 
0.08(stat.)±0.03(syst.). 
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This Letter presents a study of the charge asymmetry 
in the decay B ± — » J/ifjK ± (Tr ± ), which is defined as 

A C p(B+ -f J/iJjK + (n+)) = 
N(B~ -> J/ipK-(7T-)) - N{B+ -> J/^K+(n + )) 
N(B- -> J/t/>K-(ir-)) + N(B+ ->• J/t/>K+(ir+)) ' 

A non-zero value of Acp(B + — > J/tpK + (n + )) corre- 
sponds to direct CP violation in this decay. In the 
b — > sec transition (charge conjugate states are assumed 
throughout), the tree- level and b — > s penguin amplitudes 
have a small relative weak phase, arg[— V^ s V^/Vt s V^]. 
Therefore, the standard model (SM) predicts a small 
A C p(B+ -> J/ijjK+) ~ 0.003 Q. Thus, the measure- 
ment of Acp(i? + — ► J/tjjK + ) is an important way of 
constraining those new physics models which predict an 
enhanced value of this asymmetry, up to 0.01 or higher 
l]. Most cited are the models with an extra £7(1)' gauge 
boson responsible for the flavor-changing coupling be- 
tween b and s quarks [2] and the Two-Higgs Doublet 
Model (2HDM), which introduces an extra coupling to 
the charged Higgs boson |3j. 

In b — > dec transitions, on the contrary, the relative 
phase between the tree-level and b —> d penguin dia- 
gram, &rg[—V c dV* b /VtdV t * b ], is expected to be significant 
so that direct CP-violation may be of the order of one 
percent 0, HI ■ Decays governed by b — » dec transition 
have already been explored. Recently the Belle collabo- 
ration reported large direct CP violation in B° — > D + D~ 
decays, A D+D - = +0.91 ± 0.23 ± 0.06 @], much in ex- 
cess of the SM expectation. However, this result was not 
confirmed by the BaBar collaboration, which measured 
C D+D - = -A D+D - = +0.11 ± 0.22 ± 0.07 0. Here, we 
report a complementary measurement of the direct CP 
violation asymmetry in the b — > dec transition using the 
decay B + — > J/%jjn + . 

The DO detector is described in detail elsewhere Q. 
The detector components most important for this anal- 
ysis are the central tracking and muon system. The DO 
central tracking system consists of a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT), both 
located within a 2 T superconducting solenoidal magnet. 
The muon system is located outside the calorimeters and 
consists of a layer of tracking detectors and scintillation 
trigger counters in front of 1.8 T iron toroids, followed by 
two similar layers behind the toroids [9j . The polarities of 
the solenoid and toroid are reversed regularly during data 
taking, so that the four solenoid-toroid polarity combina- 
tions are exposed to approximately the same integrated 
luminosity. The reversal of magnet polarities is essential 
to reduce the detector-related systematics in asymmetry 
measurements and is fully exploited in this study. 

The decay chain B+ -> J/if)K + (n + ) with J/tp -> 
is selected for this analysis from 2.8 fb _1 recorded 



by DO. Each muon is required to be identified by the 
muon system, to have an associated track in the central 
tracking system with at least two measurements in the 
SMT, and a transverse momentum pt^ > 1.5 GeV/c with 
respect to the beam axis. At least one of the two muons 
is required to have matching track segments both inside 
and outside the toroidal magnet. The di-muon system 
must have a reconstructed invariant mass between 2.80 
GeV/c 2 and 3.35 GeV/c 2 . An additional charged particle 
withpT > 0.5 GeV/c, total momentum above 0.7 GeV/c, 
and at least two measurements in the SMT, is selected. 
This particle is assigned the kaon mass and is required 
to have a common vertex with the two muons, with the 
X 2 of the vertex fit being less than 16 for three degrees of 
freedom. The displacement of this vertex from the pri- 
mary interaction point is required to exceed three stan- 
dard deviations in the plane perpendicular to the beam 
direction. The primary vertex of the pp interaction is 
determined for each event using the method described 
in Ref. [lfj. The average position of the beam-collision 
point is included as a constraint. 

From each set of three particles fulfilling these require- 
ments, a B + candidate is constructed. The momenta of 
the muons are corrected using the J/ip mass constraint. 
To further improve the B + selection, a likelihood ratio 
method [11| is applied. The variables chosen for this 
analysis include the lower transverse momentum of the 
two muons, the \ 2 of the B + decay vertex fit, the B + 
decay length divided by its uncertainty, the significance 
Sb of the B + track impact parameter, the transverse mo- 
mentum of the kaon, and the significance Sk of the kaon 
track impact parameter. For any track i, the significance 
is defined as Si = \J[er / p{£t)] 2 + [cl/c^l)] 2 , where er 
(ex) is the projection of the track impact parameter on 
the plane perpendicular to the beam direction (along the 
beam direction), and a(ex) is its uncertainty. The 

track of each B + is fitted assuming that it passes through 
the reconstructed vertex and is directed along the recon- 
structed B + momentum. Finally, the mass of the re- 
constructed B + candidate is constrained to the window 
4.98 < m(J/ipK) < 5.76 GeV/c 2 . 

The resulting invariant mass distribution of the J/tpK 
system is shown in Fig. [1] with the result of an un- 
binned likelihood fit to the sum of contributions from 
B — ► J/tpK, B — > J/ipiT, and B — > J/ipK* decays, as 
well as combinatorial background (BKG). The mass dis- 
tribution of the J/ipK system from the B — > J/i^K hy- 
pothesis is parameterized by a Gaussian function with 
the width depending on the momentum of the K candi- 
date. The mass distribution of the J/ipir system from the 
B — » J/ipn hypothesis is parameterized by a Gaussian 
function with the same width. It is then transformed 
into the distribution of the J /ipK system by assigning 
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FIG. 1: The J/ipK invariant mass distribution together with 
the result from the unbinned likelihood fit (the unsplitted 
sample) . 



the kaon mass to the pion. The decay B — > J/tpK* with 
— ► i^7r, where the pion is not reconstructed, pro- 
duces a broad J/ ^if mass distribution with the threshold 
near m(B) — m(ir). It is parameterized using the Monte 
Carlo simulation. The combinatorial background is de- 
scribed by an exponential function. The J /tpK, J/i/jtt, 
and J/4>K* contributions depend on the kaon momen- 
tum. The Monte Carlo simulation shows that this de- 
pendence can be modeled by the same polynomial func- 
tion with different scaling factors for J/ipK, J/tpir, and 
J/ipK* signals. The coefficients of the polynomial are 
determined from the fit. The B — > J/ijjK signal contains 
40,222 ± 242(stat.) events, while the B -> J/ipir signal 
contains 1,578± 119(stat.) events. 

To measure the charge asymmetry A between the 
J/ipK~(n~) and J/tfjK + (-K + ) final states, both physics 
and detector effects contributing to the possible imbal- 
ance of events with positive and negative kaons must be 
taken into account. One physics source of asymmetry is 
direct CP violation in the B + — > J/tpK + (ir + ) decay. In 
addition, forward-backward charge asymmetry of events 
produced in the proton-antiproton collisions can also be 
present. Detector effects can give rise to an artificial 
asymmetry if, for example, the reconstruction efficiencies 
of positive and negative particles are different. However, 
a positive particle produces the same track as a negative 
particle in the detector with reversed magnet polarity. 
Therefore, essentially all detector effects can be canceled 
by regularly reversing the magnet polarity. 

Following the method applied in Ref. [Il,[l3j], the event 
sample of Fig. Q] is divided into eight subsamples cor- 
responding to all possible combinations of the solenoid 
polarity = ±1, the sign of the pseudorapidity of the 
J/ipK system 7 = ±1, and the sign of the kaon candi- 
date charge q = ±1. In each subsample, the number n£ 7 
of the events in the contributing channels, J/ipK, J/ipir 



and J/tpK* , is obtained from the unbinned likelihood fit 
to the mass distribution m(J/ipK) using the same like- 
lihood function as for the whole sample. All parameters 
of the fits apart from the fractions of the J/ipK signal, 
the J/tpir signal, and the J/ipK* signal, are fixed to the 
values determined from the fit to the whole sample. 

The number of events in the J/ipK and J/ipir chan- 
nels for each f3^q subsample are used to disentangle the 
physics asymmetries and the detector effects. The n^ 7 
can be expressed through the physics and the detector 
asymmetries as follows [12j: 

n fh = ±Ne f, (l + qA){l + <r f A fb )(l + >yA det ) 

x (l + qP 1 A qf)j )(l + qPA q p)(l+j3'yA M ). (1) 

Here N is the total number of signal events; e 13 is the 
fraction of integrated luminosity with solenoid polarity 
(3 (e+ + e~ = 1); A is the charge asymmetry to be mea- 
sured; Afb accounts for possible forward-backward asym- 
metric B meson production; A det is the detector asym- 
metry for kaons emitted in the forward and backward 
direction; A g 0^ accounts for the change in acceptance of 
kaons of different sign bent by the solenoid in different di- 
rections; A q p is the detector asymmetry, which accounts 
for the change in the kaon reconstruction efficiency when 
the solenoid polarity is reversed; A^ accounts for any 
detector-related forward-backward asymmetries that re- 
main after the solenoid polarity flip. We apply a \ 2 fit 
of Eq. [1] to the number of events in all subsamples and 
extract all asymmetries and the total number of events in 
the J/ipK and J/i/jtt channels together with the fraction 
of events with positive solenoid polarity e + , which is con- 
strained to be the same for both channels. Results are 
presented in Table |U The charge asymmetry between 
B~ — » J/i/jK^ and B + — > J/i{iK + is measured to be 
A(J/ipK) = -0.0070 ± 0.0060, and the charge asymme- 
try between B~ — » J/ipir~ and B + — > J/-07T+ is found 
to be A(J/ipn) = -0.09 ± 0.08. The detector asymme- 
tries are all consistent with zero, since the acceptance of 
the charged particles of different sign inside the solenoid 
is the same. However, we measure these asymmetries 
directly and do not rely on assumptions. The forward- 
backward asymmetry is also consistent with zero, as ex- 
pected in the SM. 

In addition to the detector effects, the charge asym- 
metry A[B — > J/ifjK) is affected by the difference in the 
interaction cross-section of K + and K~ with the detector 
material which is due to the fact that the reaction 
K~N — > Yir (where Y are hyperons A, S etc.) has no 
K + N analog. The difference in the interaction cross sec- 
tion results in a lower reconstruction efficiency of K~ 
and a visible kaon charge asymmetry Ak between K~ 
and K + candidates, which shifts the A(J/ipK) asymme- 
try. The kaon asymmetry is measured directly in data 
by comparing the exclusive decay c — > D* + — > D°ir + , 
D° — > [i+v^K^ and its charge conjugate. It is expected 
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TABLE I: Physics and detector asymmetries for J/ipK and 
J/ipn channels. e + is constrained to be the same for both 
channels. 





J/tPK 


J/lpTT 


N 


40, 217 ±243 


1,577 ±118 


e+ 


0.5060 ± 0.0030 




A 


-0.0070 ±0.0060 


-0.09 ±0.08 


A fb 


0.0013 ±0.0060 


0.04 ±0.09 


Adet 


-0.0033 ±0.0060 


0.21 ±0.08 


Aqff-y 


-0.0050 ±0.0060 


-0.02 ±0.09 




0.0001 ±0.0060 


-0.19±0.08 


A M 


-0.0030 ±0.0060 


0.05 ±0.08 
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FIG. 2: The Am distribution (solid line) for 1.6 < m(fiK) < 
1.7 GeV/c 2 . The background distribution of events with 
wrong charge correlation of a muon and a kaon (dashed line) 
is also shown. The background normalization is required to 
give the same number of events with right and wrong charge 
correlation for 0.19 < Am < 0.22 GeV/c 2 . 



from theory that there is no CP violation in the semilep- 
tonic D° decays [HI . The possible CP-violating effects in 
B -» D*±X decays are estimated to give a negligible con- 
tribution. Therefore, the observed asymmetry is only due 
to kaon reconstruction. The decay of D* produces a clear 
peak in the mass difference, Am — m([iKir) — m{fJ,K). 
Its width depends on the mass m(nK). An example of 
the Am distribution for 1.6 < m(^iK) < 1.7 GeV/c 2 
is shown in Fig. [2l The combinatorial background un- 
der the peak is determined using events where all three 
particles (muon, kaon, and pion) have the same charge, 
and its normalization is obtained using events with large 
values of Am outside the D* peak. The number of 
D* — > D°tt decays is determined by subtracting the nor- 
malized number of background events from the number 
of signal events in the mass band corresponding to the 
D* peak. The width of this band is varied depending 
on the mass of the \iK system to ensure maximal signal 
significance. 

The detector charge asymmetries are disentangled 
from the kaon asymmetry using the same detector model 



J/^K+) 



FIG. 3: Dependence of the kaon asymmetry on the kaon mo- 
mentum pk in eight bins of approximately equal statistics. 
Errors are statistical. 



of Eq. [TJ To account for the momentum dependence of 
the kaon cross-section [141 ], the kaon asymmetry is mea- 
sured in different bins of kaon momentum px, as shown in 
Fig. [3] The obtained asymmetry is convoluted with the 
kaon momentum distribution in the B — > J/ipK decay 
and the resulting kaon asymmetry in the B — > J /tpK de- 
cay is found to be A K = -0.0145 ± 0.0010. Taking into 
account this value, we obtain Acp{B + - 
A{J/i)K) -A K = ±0.0075 ± 0.0061(stat.) 

The systematic uncertainty of Acp(B + — > J/i/jK^) is 
estimated as follows. The systematic uncertainty from 
the unbinned fit of the invariant mass distribution of the 
J/ipK system is estimated by varying the parameters 
fixed during the unbinned fitting in the ftjq subsamples 
by ±lcr, and is found to be 0.0002. The systematic un- 
certainty from the choice of the fitting range is found 
to be 0.0004. The shape of the J/ipK* contribution to 
the likelihood function is parameterized using the Monte 
Carlo simulation, and therefore produces an uncertainty 
in the number of signal events. We repeat the fit with dif- 
ferent models of J/ipK* contribution, including a model 
without any such contribution. The maximal deviation 
in the resulting asymmetry is found to be 0.0025, which 
is taken as the systematic uncertainty from this source. 
To measure the kaon asymmetry in the detector, we sub- 
tract the combinatorial background under the D* peak 
(see Fig. [2j dashed line). To estimate the uncertainty 
from the background definition, we required the muon 
and the pion to have different charges and repeated the 
measurement of the kaon asymmetry. The resulting de- 
viation in A C p(B+ J/tpK+) is found to be 0.0008. 
Also, the sample used to measure the kaon asymmetry 
contains a contribution of D° semileptonic decays with- 
out a charged kaon in the final state. They are taken 
into account assuming the same selection efficiency as 
the dominant D° — > [ivK decay. To find the impact of 
this assumption on the final result, we repeated the mea- 
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surement of the kaon asymmetry assuming a zero recon- 
struction efficiency for additional D° decay modes. The 
resulting deviation in A CP (B+ -> J/ipK+) is 0.0005. 
Combining all contributions in quadrature, we estimate 
the total systematic uncertainty on Acp(B + — > J/ipK + ) 
to be 0.0027, which is dominated by the uncertainty from 
the mass distribution model. 

The systematic uncertainty of Aqp{B + — * J/i{nr + ) is 
estimated similarly to that of Acp(B + — » J/ipK + ). The 
only sizable contributions are: 0.01 from the variation of 
the fitting range, and 0.02 from the mass model. The 
total systematic uncertainty is 0.03. 

In conclusion, the direct CP violating asymmetry in 
the B + — > J/ipK + decay is measured to be Acp{B + — > 
J/ipK+) = +0.0075 ±0.0061(stat.)±0.0027(syst.), which 
is consistent with the world average, Aqp(B + — > 
J/ipK+) = +0.015 ± 0.017 14], but has a factor of two 
improvement in precision, thus providing the most strin- 
gent bounds for new models predicting large values of 
A C p(B+ -> J/ipK + ). The direct CP violating asym- 
metry in the _B + — > J/ipTr + decay is measured to be 
A C p(B + -> J/07T+) = -0.09 ± 0.08(stat.)±0.03(syst.). 
Our result agrees with the previous measurements of this 
asymmetry |l4j and has a competitive precision. 
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